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The light-like cusp anomalous dimension is a universal function that controls infrared divergences 
in quite general quantum field theories. In the maximally supersymmetric Yang-Mills theory this 
function is fixed fully by integrability to the three-loop order. At four loops a non-planar correction 
appears which we obtain for the first time from a numerical computation of the Sudakov form 
factor. Key ingredients are widely applicable methods to control the number-theoretic aspects of 
the appearing integrals. Our result shows explicitly that quadratic Casimir scaling breaks down at 


four loops. 


PACS numbers: 11.25.Db, 12.38.Bx 
INTRODUCTION 


Quantum field theory is used throughout physics due 
to its unrivalled power to perform systematic compu- 
tations. Connecting the calculated quantities to ob- 
served experimental variables, however, can be compli- 
cated. Perturbation theory for the theory of the strong 
nuclear force, quantum chromodynamics (QCD), for in- 
stance can famously only be trusted on very small length 
scales al Bl. The physical detectors, however, involve 
vastly larger length scales where the perturbation the- 
ory is guaranteed to break down. A central role in the 
analysis is played by universal functions, a prime exam- 
ple of which is known as the light-like cusp anomalous 
dimension (CAD). 


The universality of the CAD makes it appear in many 
contexts, e.g. (1) By definition, it refers to the anomalous 
dimension of a Wilson loop with a light-like cusp 3, [ay 
(2) It provides the leading infrared (IR) behavior of on- 
shell amplitudes and is an essential ingredient in con- 
structing amplitudes, see e.g. Bl. (3) It determines the 
logarithmic growth of anomalous dimensions of high-spin 
Wilson operators. Such operators naturally appear in 
the operator product expansion description of deep in- 
elastic scattering processes in QCD in the so-called large 
Bjorken limit, see e.g. ld. (4) It is related to the gluon 
Regge trajectory, and can be obtained from the so-called 
BFKL equation, see e.g. i. (5) Via the AdS/CFT cor- 
respondence |8], cusped Wilson operators in the planar 
limit admit a dual description in terms of spinning strings 
in a curved background [9]. There is therefore ample di- 
rect physical motivation to compute this function. 


In N = 4 super Yang-Mills (SYM) theory, the planar 
cusp anomalous dimension is known non-perturbatively 
via the so-called BES equation based on ideas of the 
AdS/CFT correspondence and integrability mi Con- 
sistency with field-theory computations through to four 
loops has been obtained at the planar level 13, and 
also confirmed by the strong coupling two-loop string 
computations . There is an intriguing connection be- 
tween N = 4 SYM and QCD results, via the maximal 
transcendentality property: the N = 4 anomalous di- 
mensions correspond to the “leading-transcendentality” 
contribution in QCD (15) [16}. Inspired by this, the three- 
loop N = 4 results were first obtained in id using the 
three-loop QCD results (17. 

At four loops in gauge theories with matter fields in 
the adjoint representation the first correction beyond the 
planar limit enters into the CAD. This Letter contains to 
our knowledge the first computation of this quantity. In 
perturbation theory four-loop non-planar integrals have 
to date proven to be largely prohibitively complicated. 
Non-planar corrections are also hard to analyse directly 
within the AdS/CFT correspondence as they correspond 
to string loop corrections, and the role of integrability be- 
yond the planar limit is unclear. Moreover, a quite gen- 
eral conjecture has been posed by extrapolating three- 
loop results that the non-planar part of the CAD van- 
ishes in any theory. This is usually called quadratic 
Casimir scaling of the CAD, and plays an important role 
in the IR factorisation in gauge theories fig ba, This 
scaling is known to break down in planar N = 4 at strong 
coupling, as well as through instanton effects [24]; in this 
Letter the approach will be strictly perturbative. 

The minimal observable that contains the cusp anoma- 
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lous dimension is the Sudakov form factor. In maximally 
supersymmetric Yang-Mills theory this is a correlator of a 
member of the stress-tensor multiplet with two on-shell 
massless states. The first computation of the two-loop 
correction to the Sudakov form factor in N = 4 SYM 
appeared in [25]. The three-loop correction to the QCD 
result was studied in a series of papers [26 {30}. In 
these results were fine-tuned for the form factor in M = 4 
SYM to the three-loop order. The integrand for the four- 
loop Sudakov form factor in M = 4 SYM was derived 
in based on the duality between color and kinemat- 
ics, and its reduction to master integrals was presented 
in B3. Recently progress for planar corrections at four 
loops in QCD was reported in [34H36]. For the five-loop 
integrand in M = 4 SYM see [37]. 


REVIEW 
Form factor and cusp anomalous dimension 


The Sudakov form factor involves only a single scale q? 
which is the Lorentzian norm of the sum of the two on- 
shell state momenta, i.e. g? = (pı + p2)? with p? = p3 = 
0. Dimensional analysis and maximal supersymmetry fix 
the form factor F™ at l loops to be given by 


FO = Peg (ag?) FO (1) 


where the coupling constant is normalised as g? = 
Run (Are”iE)f. 
(UV) divergences since the operator is protected, leaving 
only IR divergences. If dimensional regularisation with 
D = 4 — 2% is used to regulate the latter, F is a purely 
numerical function of gauge group invariants and e. This 


The form factor has no ultraviolet 


function is related to the cusp anomalous dimension Ye 
at L loops by [38 : 


l I 
] F A Ate Goi Fi (1) O 2 
(log F)“? = (le)? + De +Fin”’|+0(ed). (2) 


At I loops the planar (i.e. N! leading-color in SU (N.) 
gauge theory) part of F has leading divergence œ 1/e?. 
This function needs to be expanded down to e7? to ex- 
tract the l-loop CAD, and also higher terms in the Lau- 
rent expansion in e from lower-loop contributions are re- 
quired. The gauge group dependence of the CAD shows 
that the first occurrence of non-planar (i.e. subleading- 
color) corrections is at four loops, see e.g. , due to the 
appearance of a quartic Casimir invariant. This invari- 
ant therefore breaks quadratic Casimir scaling explicitly. 
The relation between form factor and cusp anomalous 
dimension for the non-planar part at four loops is 


(4) 


= _ Yeusp, NP As oO (e71) , (3) 


[Fo NP (8e)? 


i.e. [F®]np has only a double pole in e. However, indi- 
vidual integrals that contribute to [F“ np will typically 
show the full 1/e® divergence. 

The general CAD is expressible as a rational-coefficient 
polynomial of Riemann Zeta values ¢,, or their multi- 
index generalizations, Cn, .ns,..., known as multiple zeta 
values (MZVs) (see e.g. uA). These numbers have a 
transcendentality degree which is the sum of their indices, 
X; ni. At l loops, the planar CAD in M = 4 SYM has 
uniform transcendentality degree 21 — 2. At four loops 
for instance, the planar CAD in N = 4 SYM has been 
computed ENE) to be 


oy = 175266 — 64¢2. (4) 


We will provide strong evidence that also the non-planar 
CAD is of uniform transcendentality six at four loops. 

In QCD, the CAD has the same maximal transcen- 
dentality degree as in N = 4, but also contains lower 
transcendentality degree constants. The coefficient of the 
maximal transcendentality coefficients match between 
planar N = 4 and QCD, an observation known as the 
maximal transcendentality principle. 


Integrands, integrals, integral relations 


The non-planar part of the Sudakov form factor in 
N = 4 SYM was obtained as a linear combination 
of a number of four-loop integrals in using color- 
kinematics duality [43] hd. The integrals take the generic 
form 


I= (P? f Pu sad la TE D,’ (5) 
k=1 “k 


where D; are propagators and the numerators N (l;, p;) 
are quadratic polynomials of Lorentz products of the four 
independent loop and two independent external on-shell 
momenta. The explicit form of the integrals for the prob- 
lem at hand can be found in BA. There are 14 distinct in- 
tegral topologies that contribute to the non-planar CAD, 
labelled (21) — (34) in the given reference, each with 12 
internal lines. These coordinates and the full parametri- 
sation of the integrals as in will be used throughout. 

The integrands are identified up to terms that inte- 
grate to zero. A special class of these are generated by 
infinitesimal linear reparametrisations of the loop mo- 
menta and are known as integration-by-parts (IBP) iden- 
tities [46]. Using these identities the form factor was 
simplified in using the Reduze code (47. A particular 
subset of such relations, dubbed ‘rational IBP’ relations 
and obtained in ad, will play an important role for the 
problem at hand. Note that integral relations due to 
graph symmetries are a particular subset of the rational 
IBP relations. Although simpler integrals emerged in 
compared to BA, these have evaded integration so far due 
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to their overwhelming complexity. The first obstacle to 
computing the CAD is therefore to find a complete set of 
integrals which are each simple enough to integrate. 


UNIFORMLY TRANSCENDENTAL INTEGRALS 


The key idea is to find a representation of the four-loop 
form factor such that all integrals have uniform transcen- 
dentality if the dimensional regularisation parameter is 
assigned transcendentality —1. Such integrals will be re- 
ferred to as UT integrals. Such a representation of the 
form factor makes manifest the expected maximal tran- 
scendentality property of planar N = 4 SYM, and has 
been achieved at three loops in 31]. Crucially for our 
purposes, the UT integrals turn out to be much simpler 
to integrate than generic integrals. 

Beyond explicit computation, there are various ways 
to prove a general integral is a UT integral: 


e A UT integral can be expressed in the so-called 
dLog form m Bg. 


e A set of integral basis that are all UT integrals can 
lead to certain simple differential equations 51]. 


e The leading singularities, or equivalently, the 
residues at all poles of a UT integral must always 
be a constant Bd. 


We have obtained a dLog form for instance in topolo- 
gies (21) and (23) for the integrals with numerator 


[is - p1)?]? (6) 


each. The dLog form of the topology-(21) integral as well 
as its analytic result were also obtained in . Another 
example where it was possible to derive a dLog form is 
topology (28) with numerator 


(ls — l4 — p2)” x (l3 - pı)”. (7) 


Finding a dLog form in general, however, is a difficult 
problem. Moreover, the differential equations method 
is not directly applicable for the Sudakov form factor 
since it is a single-scale problem (see 34,53 for a work- 
around). The last criterion in the above bullet-point 
list, however, allows an algorithm. For this the four- 
dimensional loop momenta are parametrised such that 
only scalar integration parameters remain, e.g. as 


lı = aıpı + @2p2 + a3q1 + auge , (8) 


where pı and pa are the external on-shell momenta, and 
qı, q2 can be chosen as the two complex solutions to 


a =G Pj = 0 Vi, j and di" @2 = Pi “P2 (9) 


With this, the integration volume changes to 


fen en orm) [Th (...). (10) 


In this so-called parametric form, the leading singular- 
ity is studied by computing subsequent residues for all 
occurring scalar parameters (16 at four loops). For this, 
one needs to choose an order of taking residues in the 
scalar parameters. If, after taking a residue in a partic- 
ular parameter, one encounters other than a simple pole 
in a remaining parameter, the integral is not UT. In prin- 
ciple, there are 16! ~ 2 x 101° different orders in which 
the residues can be taken. We will explain below how we 
deal with this in practice. 

Besides checking that an integral at hand is UT or 
not, this procedure can be used to constrain the space of 
potential UT integrals of a given topology. We can start 
with an Ansatz of a linear combination of momentum 
products of mass dimension four. Given a set of integrals 
for a given topology one can see if one or more members 
fail the UT test for a particular randomly chosen order 
of residues. If a quadratic or higher pole is found, one 
can derive a linear constraint on the set of coefficients in 
the Ansatz to evade the pole. Solving these constraints 
gives a smaller set of integrals. In the case at hand about 
a thousand such random checks typically yields a list of 
candidate UT integrals that does not easily yield further 
constraints. Note a related strategy was used in 50). 
Through this strategy a set of candidate integrals was 
obtained for each topology. 

For ease of integration it is advantageous to have a 
form of the numerator as a product of two factors, both 
quadratic in momenta. To find linear combinations of 
UT candidate integrals with this property one writes 
a product-type Ansatz and derives a set of quadratic- 
constraint equations. In some cases, a solution may not 
exist or may not be useful. For the non-planar four-loop 
form factor, this happens in topologies (26) and (25), re- 
spectively. Here the Ansatz can be widened by adding 
integrals with two fewer propagators (i.e. two pinched 
lines). 

An example result is the topology (25) integral with 
numerator 


[(p1 — 15)? + 2(l4 — Is)? + (Ig — la)? — (lz — 15)? 
(pı la)” 4 (14 — ts)? (pı — l3 +l4—15)?. (11) 


This integral captures the full 12-line context of the non- 
planar form factor in this topology. 

Having obtained integrals which are product-form UT 
candidates — possibly up to a ten-line contribution — one 
needs to express the form factor as a linear combina- 
tion of them. Since the form factor is expected to be 
uniformly transcendental itself, this linear combination 
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must involve only rational numbers — not e. The afore- 
mentioned ‘rational IBP’ relations obtained in prove 
extremely useful for this task. 

The first result of this article is that the non-planar 
part of the four-loop form factor obtained in can be 
expressed in terms of 23 UT candidate integrals only, 
providing supporting evidence for maximal transcenden- 
tality in the non-planar sector. 

The linear combination of product-form UT integrals 
that builds the non-planar four-loop form factor is not 
unique owing to the rational IBP relations between in- 
tegrals. Our choice of basis focuses on minimising the 
number of different product-type integrals, and is guided 
by the condition that the appearing integrals can be com- 
puted numerically. The full list of all used UT candidate 
integrals including their results will be provided in an 
accompanying long write-up 53. 

Note that besides those integrals for which an explicit 
dLog form was obtained, all UT candidates have passed 
a significant number of residue checks, e.g. the two most 
complicated ones are in topology (26) and have passed 
10° such checks. This provides strong evidence that they 
are indeed UT integrals. 


NUMERICAL INTEGRATION 


There is a variety of techniques available to com- 
pute loop integrals both analytically and numerically, see 
e.g. 54) for an overview. Here mainly sector decomposi- 
tion and Mellin-Barnes (MB) integration have been used, 
as well as an exact integral result of topology (21) Bd. 

Two computer implementations of sector decomposi- 
tion have been used: mostly FIESTA bad. wita 
cross-checks for simpler integrals using SecDec |60 i 
The numerical integration is performed using the VE- 
GAS algorithm [63] as implemented in the CUBA library 

. To check on the reported uncertainties, several runs 
with different integration settings were performed. We 
observe that error bars reported by FIESTA scale with 
1/VN,, with N, the number of sampling points (typically 
in the range of a few up to several 100 millions). More- 
over, uncertainties given by FIESTA tend to be rather 
over- than underestimated which can be seen by com- 
paring to MB or analytic results, and by the fact that 
fluctuations upon increasing N, are well within the re- 
ported error bars. This error will therefore be interpreted 
as the standard deviation of a Gaußian distribution, and 
errors are added in quadrature. Further details will be 
provided in [53]. A feature of the numerical integration 
that can be observed is that UT integrals usually gener- 
ate considerably fewer integration terms within FIESTA 
compared to non-UT siblings of comparable complexity, 
even though for instance single-sector contributions are 
not UT separately. 

A second integration technique applied here is MB in- 


tegration. Several automated tools exist, such as MB 
tools [63) led and AMBRE leHeg. Given an integral 
of the form (5), one can derive a wide variety of MB 
representations, and the main challenge is to derive valid 
MB representations for crossed four-loop topologies. But 
even if the latter is possible, these representations often 
feature prohibitively complicated integrals for the prob- 
lem at hand. Here efficient MB representations for a sub- 
set of (planar and crossed) integrals have been derived by 
means of a hybrid of the loop-by-loop approach and using 
the F and U graph polynomials (we refer again to 
for details). If an efficient MB representation is avail- 
able, the precision of the integrals is typically three to 
four orders of magnitude better than FIESTA. We have 
cross-checked all our MB results with FIESTA. 


Results 


Our results for the form factor up to an overall fac- 
tor are summarised in table [I] As expected on physical 
grounds, the first six orders of expansion of the non- 
planar form factor cancel well within error bars. 


TABLE I. Non-planar form factor. 


result 


For the leading non-trivial poles, e18:6,5,4} the preci- 
sion in each integral is good enough to write the number 
as a rational multiple of {1, (a, ¢3, C4}, respectively. Note 
the e`” coefficient should vanish for each integral, and 
indeed does well within error bars. The coefficients of 
e7 {8:6:543} can be non-zero in individual integrals but 
must cancel in the linear combination of the form factor, 
which is indeed the case. We emphasize that the vanish- 
ing of all these poles provides a very strong consistency 
check for the computation, and gives further support for 
our error analysis. 

The non-trivial coefficient for the non-planar form fac- 
tor at order e7? yields the sought-after non-planar four- 
loop CAD: 


np = —3072 x (1.5931 + 0.3362) 7 (12) 
where 3072 = 2 x 24 x 64 is the normalization factor in- 
cluding that from the permutational sum, the color fac- 
tor Bd, and the denominator of (8), respectively. This 
is the second major result of this paper. The naive sig- 
nificance of a deviation from zero is 4.70. Note that 
even with a more conservative error estimate obtained 
by adding uncertainties from the individual UT integrals 
linearly, we still have 1.5931+0.9813 which is significantly 
non-zero. 
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Compared to the planar result (4): An ~ —1875, 
we see that the non-planar CAD has the same sign. If 
N. = 3 is used, its central value becomes 7s, Np ~ 
—544, i.e. the planar contribution is a factor of 3 — 4 


larger. 


DISCUSSION 


In this Letter we present the first computation of 
the non-planar correction to the cusp anomalous dimen- 
sion in the maximally supersymmetric Yang-Mills theory, 
which first starts at four loops. 

While integrating a generic set of four-loop form fac- 
tor master integral remains quite challenging, a key idea 
leading to the present breakthrough is to express the 
form factor in terms of uniformly transcendental inte- 
grals. The UT basis can be constructed algorithmically. 
While a generic integral contains mixed transcendentality 
degrees, the UT integrals contain only numbers of fixed 
(maximal for N = 4 SYM) transcendentality in each or- 
der of the e expansion. What is interesting and deserves 
further study is that such simplicity is inherited in sector 
decomposition, where each sector integral is no longer UT 
separately. Once a good set of (candidate) UT basis inte- 
grals is determined, the full form factor can be expressed 
in this basis by using a simple subset of the IBP relations 
or, in principle, directly by unitarity cuts. In the case at 
hand, our results provide strong evidence for the max- 
imal transcendentality of the four-loop non-planar form 
factor in MN = 4 SYM. In particular this implies strongly 
that the four-loop CAD can be written as a rational lin- 
ear sum of (g and (2. The UT basis finding algorithm is 
widely applicable. An immediate interesting application 
of already obtained results would be to four-loop propa- 
gator integrals in QCD, see e.g. [zo]. A UT basis in QCD 
would always involve pre-factors with non-negative pow- 
ers of e, simplifying computations potentially drastically, 
see for similar ideas. 

Based on the UT basis, a numerical computation has 
also yielded the first information on the value of the CAD 
in N = 4, which is statistically significantly non-zero. In 
particular, our result indicates that quadratic Casimir 
scaling breaks down at four loops in this theory and is 
therefore not expected to hold in any other theory such as 
QCD. Direct further research directions include improv- 
ing precision at order 1/e? and computing further orders 
in the «expansion which contain the so-called collinear 
anomalous dimension. Computing the non-planar CAD 
in other theories such as QCD at four loops is a prime 
further target. It would be extraordinarily interesting to 
extend integrability arguments to the non-planar CAD 
in M = 4; our weak-coupling result provides the first 
data-point for this effort. 

A more general lesson is that by focusing on the 
number-theoretically simple integrals we obtained a set of 


these that can be integrated relatively easily using meth- 
ods that at face value have little to do with number the- 
ory. In other words, simplicity is persistent in Feynman 
integrals. 
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